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The fastest possible continued fraction approximations of a class 

of functions 

Xiaodong Cao*, Yoshio Tanigawa and Wenguang Zhai 


Abstract 

The goal of this paper is to formulate a systematical method for constructing the fastest 
possible continued fraction approximations of a class of functions. The main tools are the 
multiple-correction method, the generalized Mortici’s lemma and the Mortici-transformation. 
As applications, we will present some sharp inequalities, and the continued fraction expan¬ 
sions associated to the volume of the unit ball. In addition, we obtain a new continued 
fraction expansion of Ramanujan for a ratio of the gamma functions, which is showed to be 
the fastest possible. Finally, three conjectures are proposed. 


1 Introduction 

Let f{x) be a function defined on (0, -boo) to be approximated. We suppose that there exists a 
fixed positive integer v and a constant c 7 ^ 0 such that 

( 1 . 1 ) lim x'^fix) = c. 

^ a;->-+oo ^ 

* Corresponding author. 

E-mail address', caoxiaodong@bipt.edu.cn (X.D. Cao), tanigawa@math.nagoya-u.ac.jp (Y. Yanigawa), 
zhaiwg@hotmail.com (W.G. Zhai) 

2010 Mathematics Subject Classification :11Y60 41A25 34E05 26D15 

Key words and phrases: Continued fraction; Gamma function; Rate of convergence; Inequality; Multiple- 
correction; Volume of the unit ball. 

This work is supported by the National Natural Science Foundation of China (Grant No.11171344) and the 
Natural Science Foundation of Beijing (Grant No.1112010). 

Xiaodong Cao: Department of Mathematics and Physics, Beijing Institute of Petro-Chemical Technology, 
Beijing, 102617, P. R. China 

Yoshio Tanigawa: Graduate School of Mathematics, Nagoya University, Nagoya, 464-8602, Japan 

Wenguang Zhai: Department of Mathematics, China University of Mining and Technology, Beijing 100083, P. R. 

China 


1 



In this case, we say that the function /(x) is of order x ^ when x tends to infinity, and 
denote 

(1.2) R(/(x)) := i., 

where v is the exponent of For convenience, R(0) is stipulated to be infinity. Hence, R(/(x)) 
characterizes the rate of convergence for /(x) as x tends to infinity. Prom (jl.ip . there exists a 
large positive number Xq such that f{x)/c > 0 when x > Xq. 

In analysis, approximation theory, applied mathematics, etc., we often need to investigate 
the rational function approximation problem. Let be an approximation to /(x) as x 

tends to inhnity, where F)(x) and Qm(x) are polynomials in x. Quite similarly to the rational 
approximation problem for an irrational number, in order to find a better approximation to /(x), 
we have to increase the degrees of both Pi{x) and Qm{x). The main interest in this paper is to 
try to look for the fastest possible continued fraction approximation or guess its approximation 
structure for /(x) as x tends to infinity. 

The paper is organized as follows. In Sec. 2, we mainly introduce a definition to classify 
the continued fraction. In Sec. 3, we will prepare two preliminary lemmas for later use. In 
Sec. 4, we first develop further the previous multiple-correction method. Secondly, we intro¬ 
duce a transformation named as Moritici-transformation to change a kind of continued fraction 
approximation problem. In addition, we also give its Mathematica program for the reader’s 
convenience. Thirdly, similarly to Taylor’s formula, we introduce two definitions of the formal 
Type-I and Type-II continued fraction approximation of order k for a function, and the formal 
continued fraction expansion, respectively. This section constitutes the main part of this paper. 
To illustrate our method formulated in Sec. 4, in Sec. 5 we use the volume of the unit ball as an 
example to present some new inequalities. In Sec. 6, we test the well-known generalized Lord 
Brouncker’s continued fraction formula, and show that it is the fastest possible. We also give 
some applications for the continued fraction formula involving the volume of the unit ball. In 
Sec. 7, we will use a continued fraction formula of Ramanujan to illustrate how to get the fastest 
possible form of the continued fraction expression. In Sec. 8, we explain how to guess the fastest 
possible continued fraction expansions, and give three conjectures associated to the special rate 
of gamma functions. In the last section, we analyze the related perspective of research in this 
direction. 

2 Notation and definition 

Throughout the paper, we use the notation [xj to denote the largest integer not exceeding x. 
The notation Pfc(x)(or Qk{x)) means a polynomial of degree k in x. We will use the <h(A;;x) 
to denote a polynomial of degree A: in x with the leading coefficient equals one, which may be 
different at each occurrence. While, the notation '!'(/?; x) means a polynomial of degree k m. x 
with all coefficients non-negative, which may be different at each occurrence. Let {an)n>i and 
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{bn)n>o be two sequences of real numbers with 0^/0 for all n G N . The generalized continued 
fraction 


( 2 . 1 ) 


T = bo + 


ai 


bi + 


Q2 


= bo 


Ol 


Q 2 

+ b2 +■ 


= bo+K ^ 


=1 V h 


62 + 


is defined as the limit of the nth approximant 


( 2 . 2 ) 


A.JJ 

Br, 


= 60 + K 




k=l \ bk 


as n tends to infinity. The canonical numerators An and denominators Bn of the approximants 
satisfy the recurrence relations (see [ 8 , p. 105]) 


(2.3) 


^n+2 — bn+2-^n+l T On+2^n) Bn+2 — bn+2Bn+l T (^n+2Bn 


with the initial values Aq = bo, Bo = 1 , = 60^1 + oi and Bi = 61 . 


To describe our method clearly, we will introduce two definitions as follows. 

Definition 1. Let cq 7 ^ 0, and a; be a free variable. Let {an)^=o, {bn)^=o and (cn)^o be three 
real sequences. The formal continued fraction 


Co 


(2-4) - .,^0, 

$(n;x)+K„=o Vi+ 6 ; 

is said to be a Type-I continued fraction. While, 


Co 


(2-5) - 

^u-,x) + tin=o\ ,,-2^bn.+c. 

is said to be a Type-II continued fraction. 


Remark 1. The Type-I and Type-II are two kinds of fundamental structures we often meet. 
Certainly, we may define other-type continued fraction. Because of their complexity, in this 
paper we will not discuss the involved problems. 


Definition 2. If the sequence {bn)^=o is a constant sequence (5 )^q in the Type-I ( or Type-II) 
continued fraction, we call the number uj = b (or to = |) the MC-point for the corresponding 
continued fraction. We use x = x -|- to to denote the MC-shift of x. 


If there exists the MC-point, we have the following simplified form 


( 2 . 6 ) 


Co 


$i(i/;x) + K“o (x) 



where = Cn - 


3 













3 Two preliminary lemmas 


Mortici [23] established a very useful tool for measuring the rate of convergence, which claims 
that a sequence {xn)n>i converging to zero is the fastest possible when the difference {xn — 
Xn+i)n>i is the fastest possible. Since then, Mortici’s lemma has been effectively applied in 
many papers such as UnmilliaillllMlESlEelET]. The following lemma is a generalization of 
Mortici’s lemma. For details, readers may refer to |12j . 

Lemma 1. If \mix^+oo f {x) = 0, and there exists the limit 

(3.1) lim (/(x) —/(x + 1)) = / e M, 

21—>-+00 


with A > 1, then 
(3.2) 


lim x^ ^fix) 
21—>-+00 


I 

A- 1' 


In this paper, we will use the following simple inequality, which is a consequence of Hermite- 
Hadamard inequality. 

Lemma 2. Let f be twice differentiable with f" continuous. If f"{x) > 0, then 


(3.3) 



f{x)dx > /(a + 1/2). 


4 The multiple-correction, the Mortici-transformation and the 
formal continued fraction expansion 

4.1 The multiple-correction method 

In this subsection, we will develop further the previous multiple-correction method formulated 
in mi [III. For some applications of this method, reader may refer to nnunKnun]. In fact, 
the multiple-correction method is a recursive algorithm, and one of its advantages is that by 
repeating correction-process we always can accelerate the convergence. More precisely, every 
non-zero coefficient plays an important role in accelerating the convergence. The multiple- 
correction method consists of the following several steps. 

(Step 1) The initial-correction. The initial-correction is vital. Determine the initial- 
correction <ho(;z;x) such that 


(4.1) 


R /(x)- 


4>o(i^;x)y 


= max R ( /(x) — 


<I)(z/; x) 
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(Step 2) The first-correction. If there exists a real number kq such that 


(4.2) 


R fix) - 


^oiiy,x) + ^ 

then we take the first-correction MCi(x) = ''^ith 


> R fix) - 


4>o(i^;x)y ’ 


(4.3) 


Ao = maxR fix) — -r-— 


In this case, the first-correction has the form Type-I. Otherwise, we take the first-correction 
MCi(x) in the form Type-If i.e. MCi(3:) = ^ such that 


(4.4) 


(Ato, Ao,i, Ao, 2 ) = max R ( fix) - ^ ^ ^- 

«,Ai,A 2 y ^QiV,X) + ^2+AiX-hA2 


If kq = 0, we stop the correction-process, which means that the rate of convergence can not be 
further improved only by making use of Type-I or Type-II continued fraction structure. 

(Step 3) The second-correction to the /cth-correction. If MCi(x) has the form Type-I, 
we take the second-correction 


(4.5) 

which satisfies 

(4.6) 


MC 2 (x) = 


KO 


X -|- Aq “1“ 


K 1 

x-l-Ai 


(ki,Ai) = ma,xR fix) - —-—- — 


a;-|-A 


Similarly to the first-correction, if k\ = 0, we stop the correction-process, 
ff MCi (x) has the form Type-II, we take the second-correction 

Kq 


(4.7) 

such that 


MC 2 (x) = 


x2 -h Ao,ix Ao,2 + 


(4.8) 


(ki, Ai,i, Ai, 2 ) = rnax R fix) - —-- 

k,Ai,A 2 \ 9oiiy; x)-\--rj 


x2+Ao iX-l-Ao 2-|-^7—^ 


If Ki = 0, we also need to stop the correction-process. 

If we can continue the above correction-process to determine the feth-correction function 
MCfc(x) until some k* you want, then one may use a recurrence relation to determine MCfc(x). 
More precisely, in the case of Type-I we choose 


(4.9) 


fc-i 

MCfc(x) = K , 

’ j=o\x + \j 
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such that 
(4.10) 


(Kfc_i, Afc_i) = maxR 

K.,X 


fix) - 


Kq 


K'k-2 


4“ X + Aq +•••+ X + Afc _2 + x + A 


While, in the case of Type-II we take 


(4.11) 

which satisfies 


k-l 

MCfc(x) = K 


1=0 V x^ + Xj ix + A 


1,2 


(4.12) 

where 


{Kk-i,Xk-i,i,Xk-i,2) = max R(/(x) - G(fi:, Ai, A 2 ;x)) , 

K,Ai,A2 


G(k, Ai, A2 ;x) 


c Kq _ Kk -2 _ K 

4>o(r';x) + x^ + Ao,ix + Ao ,2 +•••+ x^ + Afc_ 2 ,ix + Xk— 2,2 + x^ + Aix + A 2 


Note that in the case of both Type-I and Type-II continued fraction approximation, if Kk-i = 
0, we must stop the correction-process. In other words, to improve the rate of convergence, we 
need to choose some more complex continued fraction structure instead of it. 

Remark 2. Sometimes, we need to consider its equivalent forms. For example, the Stirling’s 
formula reads (See, e.g. [1, p. 253]) 

(4.13) r(x -|- 1) ~ V^TTx ^ , X —>■-|-oo, 
which is equivalent to 

(4.14) lim x^/(x) = 1, 

X^OO 

where 

(4.15) /(x) = Svr^ ^ r“®(x-|-l). 

From the above asymptotic formula, we may study Ramanujan-type continued fraction approx¬ 
imation for the gamma function. For more details, see Cao [12] or next section. Moreover, we 
note that ()4.13l) has many equivalent forms. Hence, it is not difficult to see that the equiv¬ 
alent transformation of a practical problem influences directly the initial-correction and final 
continued fraction approximation. 

Remark 3. If is a negative integer, our method is still efficient, i.e. we may consider the 
reciprocal of /(x). 

Remark 4. For comparison, we use the mathematical notation “R” and “max” in the above 
definition, which make the method more clearly. 
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4.2 The Mortici-transformation 

In this subsection we will explain how to look for all the related coefficients in $ 0 ( 1 ^; x) and 
MCfc(x). If we can expand /(x) into a power series in terms of 1/x easily, then it is not difficult 
to determine <I*o(^j x) and MCfc(x). Similarly, if we may expand the difference f{x) — f{x+l) into 
a power series in terms of 1/x, by the generalized Moritici’s lemma we also can find d>o(zz;x) 
and MCfc(x), e.g. the Euler-Mascheroni constant, the constants of Landau, the constants of 
Lebesgue, etc. (See [H]). However, in many cases the previous two approaches are not very 
efficient, e.g. gamma function (see, Remark 2) and the ratio of the gamma functions (for 
example, see Sec. 7 below). Instead, we may employ the following method to achieve it. 


First, we introduce the fcth-correction relative error sequence {Ei^{x))k>o follows 


(4.16) 

(4.17) 


fix) = 
fix) = 


4>o(z^;x) 


exp(Eo(x)), 


4>o(i^;x) +MCfc(x) 


exp(£’fe(x)). A: > 1, 


where 4o(A;;x) is a polynomial of degree v m. x with the leading coefficient equals one, to be 
specihed below. 

It is easy to verify that 


fix) - 
fix) - 


4>o(i^;x) 4>o(i^;x; 

c 


(exp (Eo(x)) - 1), 
c 


4>o(i^;x) + MCfc(x) 4>o(i^;x)+MCfc(x) 

It is well-known that 


(exp (^^(x)) - 1), A: > 1. 


by lima;_>.oo Ek (x) = 0 we obtain 


(4.18) 

(4.19) 


R /(x) - 


4>o(z^;x) 


R /(x) - 


t 


= V + R{Eo{x )), 

= zz -|- R (Efc(x)), A; > 1. 


4>o(i^;x) MCfc(x) 


In this way, we turn the problem to solve R(£'fc(x)). 

Take the logarithm of (I4.16h and (I4.17jl . respectively, we deduce that 

In = - In (4>o(z/; x)) -h E'o(x), 
c 

In = - In ($o(j^; x) + MCfc(x)) -|- Ek{x), k > 1. 
c 
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Next, let us consider the difference 


(4.20) 

(4.21) 


E„(x)-E„(x + l)=ln,^+ln. 


f{x + l) 4>o(i^;x + l)’ 

IT. ^ N IT. I ,14 1 fi^) , 1 4>o(z^;x)+MCfc(x) 

Ek{x) — Ek[x + 1) = In ——TT + In ' 


f{x + l) 4>o(i^;x +1)+MCfc(x +1)’ 
By Lemma 1 (the generalized Moritici’s lemma), we have 
(4.22) R {Ekix)) = R (Ekix) - Ek{x + 1)) - 1. 

Finally, if set MCo(a;) = 0, then we attain the following useful tool. 


k> 1. 


Lemma 3. Let f{x) satisfy (II.ip . Under the above notation, we have 

" 4>o(i';i) + MCt(i) 

=,,_l+R|ln-/hI- + ln ^ I, k>0. 

/(x + 1 ) <ho(^) ® + 1 ) + MCfc(3^ + 1 ) 


The idea of Lemma 3 is first originated from Mortici [23], which will be called a Mortici- 
transformation. We would like to stress that Mortici-transformation implies the following as¬ 
sertion 


, k > 0. 


(4.24) max R f{x) — - -r- 

k,A (or h:,Ai,A 2 ) \ ^o(^? ^) 4“ MC^(x) 

13 M ,1 4>o(z^;x)+MCfc(x) 

max R I in - tt + hr ■ 


k,A (or k,Ai,A2) \ /(^ T 1) ^o(^) X + 1) + MCfc(x + 1) 


In the seqnel, we will use this relation many times. For the sake of simplicity, we will always 
assume that the difference 


(4.25) 


/(!/=) 


c c f{l/z-\-l) 

is an analytic function in a neighborhood of point z = 0. 

For the reader’s convenience, we would like to give the complete Mathematica program for 
finding all the coefficients in $o(2^j3:^) and MCfc(x) by making nse of Mortici-transformation. 


(i). First, let the fnnction MT[x] be defined by 


MT[x] := In 


fix) 


In 


4>o(i^;x) MCfc(x) 


f{x-\-l) 4>o(i^;a; + 1 )MCfc(x + 1)' 

















(ii) . Then we manipulate the following Mathematica command to expand MT[x] into a 
power series in terms of 1 /x: 

(4.26) Normal[Series[MT[x]/. l/u,{u,0,lk}]]/■ Ijx (// Simplify) 

We remark that the variable 1^ needs to be suitable chosen according to the different function. 

(iii) . Taking out the first some coefficients in the above power series, then we enforce them 
to be zero, and finally solve the related coefficients successively. 

Remark 5. Actually, once we have found MCfc(x), (j4.26p can be used again to determine the 
rate of convergence. In addition, we can apply it to check the general term formula for MCfc(x). 


4.3 The formal continued fraction expansion 


Similarly to Taylor’s formula, if the A:th-correction MCfc(x) for f{x) has the Type-I (or the 
Type-II) structure, then we may construct the formal Type-I (or Type-II) continued fraction 
approximation of order k for f{x) as follows: 


(4.27) 


CFkifix)) 


1 

4>o(z^;a;) +MCfc(x)’ 


A; > 0. 


For example, Euler-Mascheroni constant has the formal Type-I continued fraction approx¬ 
imation of order k, while both Landau’s constants and Lebesgue’s constants have the formal 
Type-II continued fraction approximation of order k. For details, readers may refer to Ha- 


Example 1. Let f{x) = Then CFk{f{x)) is the Type-I, its MC-point w equals 

i (i.e. \ra = |), and 


(4.28) 


CFk{f{x)) = 


{x + |)^ + i|8(^ + l) + MCfc(a;) 

wtiprp (k K K 1 — f — 1483 8923253 10136617131375 2313439127848201 1 

wneie ai, ■ ■ ■) 32758 , 2688 ’ 159452 16 ’ 67 75390309888 ’ 1428763287038592 ’ • • • 

Example 2. Let Gn{x) be defined by (|8.22l) below. In the case of r/ = ^, CF^ (G^(x)) is 
the Type-II, for details see Corollary 2 in Sec. 7. If r/ / then CFk (G^(x)) is the Type-I, and 
it has not MC-point. We have 


(4.29) 


{Glix)) = 


-|- 2r/(l — nfjx + 2rf‘{r] — 1)^ -|- MC 2 (x) ’ 


where kq = Aq = 

V _ (2f?-i)^ I 1 I 3 
and Ai - -h 4 + 8 ( 2 r,-i)^ 


3 - J_ 

4 8(2r;-l)^ ’ 64 

6(2?7—l)'*Ki ■ 


2ri — 3) (2r/ -|- 1) -|- 10 -|- ^2r,-i)4 


If we rewrite GFk{f{x)) in a rational function of the form then s = /c -|- zz in the case 

of Type-I, and s = 2k + v m the case of Type-II. If we let R(/(x) — CFk{f{x))) = K, then 

(4.30) f{x) = GFk{f{x)) + 0{x~^), x^oo. 
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Let 00 = 0 or 1. A lot of computations reveal that if CFk{f{x)) is the Type-I, then K = 
2k + 2^ + 1 + 6 q , and AT = 4A: + + 1 + 0o hr the case of Type-II, respectively. 

For a suitable “not very large” positive integer k, by using of Mortici-transfomation and (I4.26p . 
we may get the rate of convergence for f{x) — CF]^{f{x)) when x tends to infinity. Moreover, by 
making use of telescoping method, Hermite-Hadamard inequality, etc, sometimes we can prove 
sharp double inequalities of f[x) — CFk{f{x)) for as smaller x as possible. We will give an 
example in Sec. 5. 

Now let k tend to oo, we get the formal Type-I (or Type-II) continued fraction expansion for 
f{x), or shortly write 

(4.31) f{x)^CF{f{x)):=CFooif{x)), x ^ oo. 

In some cases, we can test and guess further the general term of CF{f{x)). Here we need 
to apply some tools in number theory, difference equation, etc. We will show some examples in 
Sec. 7. 

For the formal continued fraction expansion, we are often concerned with the following two 
main problems. 

Problem 1. Determine the domains of convergence for the formal continued fraction expansion 
CF(f{x)). We may refer to two very nice books: L. Lorentzen and H. Waadeland|21]. and A. 
Cuyt, V.B. Petersen, B. Verdonk, H. Waadeland, W.B. Jones or some other classical books 
cited in there. 

Problem 2. Prove an identity for as the large domains as possible. That is, based on Problem 
1, to determine the intervals I such that f{x) = CF{f{x)) for all x G I. For example, with the 
help of continued fraction theory, hypergeometric series, etc., we hope at least to find a interval 
(xo,oo) C I for some xq > 0. Certainly, we may extend it to a complex domain. However, in 
this paper we will not investigate this topic. 

On one hand, to determine all the related coefficients, we often use an appropriate symbolic 
computation software, which needs a huge of computations. On the other hand, the exact 
expressions at each occurrence also takes a lot of space. Hence, in this paper we omit some 
related details for space limitation. 

Remark 6. From the above discussion, we observe that for a specific function, except a huge 
of computations, probably only such two kinds of structures can not provide “good continued 
fraction approximation”. In addition, in the theory of classical continued fraction, even if there 
is a continued fraction expansion for a given function, we often do not know whether it is the 
fastest possible or best possible. Generally speaking, for a given continued fraction, finding the 
rate of convergence for the A:th approximant is not always easy. 
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5 The volume of the unit ball 


It is well-known that the volume of the unit ball in M” is 


(5.1) 




n 

vr 2 

r(f+ !)■ 


Many authors have investigated the inequalities about the e.g. see [21 a n 0 El I ig [201 
Eg 12511231301E2] and references therein. 

Chen and Li [16] proved (a = |, 6 = i): 


(5.2) 


1 

Y^7r(n -I- a) 


(v)* 




1 

^/7r{n + b) 


(v)* 


Recently, Mortici [25, Theorem 3] showed that for every integer n > 3 in the left-hand side 
and n > 1 in the right-hand side, then we have the following Gosper-type inequalities: 


(5.3) 
where 

Now we let 

(5.4) 


/27re\ 




/27reY 


Y7r(n -|- 9{n)) \ n J ” \/TT{n + i9(?t,)) \ n J ’ 


N 1 1 31 

^ ^ 3 18n 810n2 


, 'i9(n) = 9{n) — 


139 

9720n3■ 


V{x) = 


TT 


r(a: + l)’ 


X > 0. 


Let us imagine that if ~ H{x) when x tends to infinity, then V{x) has an asymptotic 

formula of the form tt^H[x). In this sense, by Remark 2 and Remark 3, it suffices to consider 
the asymptotic formula for the gamma function. In fact, we note that both f[x) and l//(x) 
have the same fcth-correction MCk{x). 

From (I4.14p . we introduce the relative error sequence {Ek{x))k>o to be defined by 


(5.5) 

(5.6) 


f{x) :=87r3(^)'"r-Yx + l) 


exp(F;o(x)) 


d>o(x) 


f{x) ■■= 


exp(.Efc(x)) 

4>o(x) -bMCfc(x)’ 


A: > 1, 


where <I)o(x) = x^ + -|- |x -|- and 


(5.7) 


fc-i 

MCfc(x) = K 


j=o V X -I- 
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here kq = —i^,Ao = ^ 5^1 = nilso ’ ~ — - 1^^^98882 ’— stress that <ho(a:^) was claimed 
first by Ramanujan m, and some more coefficients may be founded in [T^]. By employing 
Lemma 1, (I5.15p (see below) and (I4.26p . it is not difficult to verify that 


(5.8) 

(5.9) 


lim x'^Eo{x) 

X—^OO 

lim x^Ei{x) 

X^OO 


11 . ^ 
1920 ■ 
459733 

~124185600 


:=Ci. 


The following theorem tells us how to improve the above results and obtain some sharper 
estimates for Eq{x) and Ei{x). 

Theorem 1. Let Eq{x) and Ei{x) he defined as (15.5p and ()5.6p . respectively. 

(i) For every real number x > 6 in the left-hand side and x > 12 in the right-hand side, we 
have 


(5.10) 


11 1 
1920 {x + 3)4 


< Eq{x) < 


11 

1920(x -5)4' 


(ii) For every real number x >9 in the left-hand side and x > 10 in the right-hand side, then 


(5.11) 


459733 1 ^ p . X ^ 459733 1 

124185600 (x - 2)6 ^ “124185600 {x + 2f' 


Proof. We use the idea of Theorem 2 in [33] or Theorem 1 in [TO] . Let Gk{x) = E^(x) — (x +1) 

for k>0. We will employ the telescoping method. It follows from Ek{x) = 0 that 


(5.12) 


Ek{x) = '^Gk{x + m), (A: = 0,1). 


m=0 


If g{oo) = g'{oo) = 0, it is not difficult to prove that 


(5.13) 


poo poo / poo \ 

g{x) = - j g'is)ds = J ij g”{t)dt \ ds. 


Note that the convenience MCo(x) = 0 . By (15.5p and (15.6K . we have 


(5.14) 

(5.15) 


By using 


Ek{x) = —61nr(x + 1 ) + 21n27r + 6 x(lnx — 1) + ln($o(a;) + MCfc(x)), 


Gk{x) = Ek{x) - Ek{x + 1) = 6 ( 1 - xln(l + -) ) + In 


d>o(g;) + MCfc(a;) 
$ 0 ( 3 ^ + 1) + MCfc(x + 1) 


Mathematica software, we can check that if x > 0, then 


(5.16) 

(5.17) 


G'o'(x) 

G'o'(x) 


11 ,29 _ 4/i(15;x) 

16x^ 8x® 96x®(l + x) 2T2(12; x) ’ 

11 , 29 9031 _ ^'3(13;x) 

16x7 ^ “ 800x9 - “800x9(1 + x)2T2(12;x) 
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9031 


By (15.131) . we get that when a: > 0, 

(5.18) “ 


480x® 336a:® 


Similarly, if x > we have 


(5.19) 


(5.20) 


and 
(5.21) 


G" (x) + 


< Go(x) < 


39872247 


11 


29 


+ 


480x5 336x6 44800x7’ 


459733 
369600x9 “ 4743200x16 

4^4(20;x)(x-^) + f||g^ 


28459200x16 4-5(6; x) (4-6(3; x)(x - ^) + 4-7(8; x) 


< 0 , 


G"(x) + 


459733 


39872247 1092949825573 

+ 


369600x9 4743200x16 32724285440x11 

4-sl20' xllx-+ ^688-"7275 

._ lel ^ 1125---2624 _ 

490864281600x11 4 - 5 ( 6 ; x) ('^6(3; x)(x — ^) + ’ 2 ^ 984 ? )^ 4 - 7 ( 8 ; x) 


> 0 , 


459733 13290749 1092949825573 ^ , 459733 13290749 

+ n < Gl(x) < — + 


20697600x7 113836800x8 2945185689600x9 

Now, combining (15.12^ . (I5.18p and (I5.2ip . we attain that 

00 ^ 00 ^ 


20697600x7 113836800x8 


(5.22) 0 < Eo{x) - ^ ^ 


1 


29 ^ 1 9031 ^ 

480 ^ (x + m)5 336 ^ (x + m)® ^ 44800 ^ (x + m)7’ 

m=0 m=0 ^ ' m=0 


(x > 0), 


, , 1092949825573 ^ 

( 5 - 23 ) - 2 ^ 


1 


2945185689600 ^ (x + mf 

m=0 ' 


^ , , 459733 ^ 


< 


20697600 ^ (x + m)7 113836800 ^ (x + m)8 

m=0 n ^ ^ 


13290749 1 „ . 1 

L 7::-^38 < 0 - (^ > 


m=0 


16' 


Let j >2 and x > i. By Lemma 2, we obtain 


(5.24) 


1 


(j - l)xl 


-1 


dt ^ 


/„ id ^ (x + m)J 

m=o ^ ’ 


< 


— /-x+m-i 


E 


dt 


_r. J 


m=0 


1 P 




lYi-l' 


By applying (5.22) and (5.24), under the condition x > 6 we have 

11 1 29 1 


(5.25) 


Eq{x) > 


480 4x4 336 5(x - \f 


11 


1 


> 


4 -1(7; x)(x - 6) + 2164192911 
■ 1920 (x +3)4 ' 13440x4(3 + x)4(-l + 2x)5 

11 1 


+ 


1920 (x + 3)4' 


13 
















































Similarly to (5.25), if x > 12, then 

11 1 


(5.26) 


Eq{x) < 


29 1 9031 

+ 


1 


480 4(x - 1)4 336 5x5 44300 6(x - \f 

11 4-2(9; x)(x - 12) + 12561000435989768 


<- 


■l920(x -5)4 
11 


67200(-l + x)4x5(-l + 2 x)6 


1920(x - 5)4’ 

This completes the proof of assertion (i). Finally, it is not difficult to check that if x > 9, then 


(5.27) F;i(x)>- 


459733 


1 


+ 


13290749 1 1092949825573 


1 


20697600 6(x- 1/2)6 113836800 7x7 2945185689600 8(x - 1/2)8 

459733 1 ^ 3 ( 13 ; x)(x - 9) + 67733478135399363858702201 


> - 


124185600 (x - 2)6 
459733 1 


736296422400(-2 + x)6x7(-l + 2x)8 


124185600 (x - 2)6’ 
and if X > 10, we have 


/r oo^ T. , ^ 459733 1 13290749 1 

(5.28) Ei{x) < - ^ , + 


20697600 6x6 113836800 7(x - 1/2)7 

459733 1 4-4(11; x)(x - 10) + 470994290293217661904 


< - 


20697600 6(x + 2)6 
459733 1 


2390572800x6(2 + x)6(-l + 2x)7 


124185600 (x + 2)6’ 

This will finish the proof of Theorem 1. 

Theorem 2. Assume n > 24, we have the following Ramanujan-type inequalities 


(5.29) 

If n> 20, then 

(5.30) 


1 /27re\2 

TT \ n J 


1 


<nn< 


1 / 27 re \2 exp (^72o(n-io)'-) 


n3 + n2 + | + i 


\ n J 6 


n3 + n2 + f+ ^ 


1 /27re 


1 - 


459733 


11642400(-4+n)« 


77 V 71 / 6 / 3 I I 21 I J__847_ 

^ ^ 2 ^ 30 9240n+9480 


< n*. 


< 


1 /27re 

^ J e/n3 + n2 + n i - 


847 


2 ' 30 9240n+9480 


□ 
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Proof. It follows from (15.ip . (15.51) and (15.61) that 

n 1 f 2TreY exp (lEoi^)) 1 /27rey exp(|.gi(f)) 

^ ^ " 0F U y y8$o(t) + 8MCi(f)- 

Now (5.29) follows from (5.10) and (5.31). 

We begin to prove (5.30). It is well-known that exp(t) > 1 -|- t. When n > 20, by the 
inequality of the right-hand side in (5.11), we have the following trivial estimate 


(5.32) 


exp 


iEi(f) 


< 1 . 


In addition, by the lower bound in (5.11), we get 


, , / 1 1 ,n, 459733 

5.33 exp - > 1-h-.El -) > 1-^re > 20 . 

^ ’ ^'6 ^27 6 ^2^ 11642400(-4 -Fre 6 ’ ^ ’ 


Combining (5.31), (5.32) and (5.33) completes the proof of (5.30). □ 

Following the same approach as Theorem 2, it is not difficult to prove the following Ramanujan- 
type inequalities for the gamma function. 


Corollary 1. Let x > 12. Then 

^ ( 8 ^“’ + + 55 ) ‘ “P (- 11520(^5)- ) < + 1) < ^ (jf (8-“ + 'I-' + + 35 

Remark 7. It should is noted that the method described in Theorem 1 and 2 also can be used 
to look for CFk{F{x)), and prove some inequalities involving the ratio of gamma functions. 

Remark 8. We will give some other results involving Lin in the subsection 6.2. 


6 Lord Brouncker’s continued fraction formula 

6.1 Lord Brouncker’s continued fraction formula 

The following formula is taken from Corollary 1 of Berndt [8, p. 145], which was first proved by 
Bauer [7] in 1872. 


Lemma 4. //Rex > 0, then 


r^(l(^ + i)) 


( 6 . 1 ) 


A °° 

^ K 


r 2 (i(x-/ 3 )) X+m=l 


(2m — 1)^ 
2 x 


By taking x = 4re -|- 1 in the above formula, we obtain the so-called Lord Brouncker’s 
continued fraction formula 


( 6 . 2 ) 


q{n) := 


r^(n + |) _ 

r 2 (re-kl) 4 ^ + 1 + 




2(4n+l)-r- 


2(4n+l)+- 


2(4n+l) + 
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For a very interesting history of formula (|6.ip . see Berndt [8, p. 145]. In addition, Lord 
Brouncker’s continued fraction formula also plays an important role in Landau’s constants, 

see [inilll]. 

The main aim in this subsection is to illustrate (without proof) that the formula (16.1|) is 
the fastest possible by making use of the method formulated in Sec. 4. Replacing x by 4x + 1 
in (16.11) and then making some simple calculation, we obtain its equivalent forms as follows. 

Lemma 5. Let Rex > — we have 


(6.3) 


r2(x + i) 


x + i + 


a2 


k: 


m + l)^ 
(i4 


Proof. Now, we are in a position to treat the above formula directly. Let 


(6.4) 


fix) 


r2(x + i) 
r2(x + i)' 


By the recurrence relation r(x + 1) = xr(x), we have 


(6.5) 


fix) ^ jx + 1)^ 
/(x + 1) (x + l)2' 


By the Stirling’s formula, it is not difficult to prove 


( 6 . 6 ) 


b_ar(x + o) 
hm x" --4 

x^oo r(x + o) 


= 1 . 


Also see [1, p. 257, Eq. 6.1.47] or [8, p. 71, Lemma 2]. 
It follows readily from (16.61) that 


(6.7) lim x/(x) = 1, 

21—>-+00 

i.e, we take v = 1 ii\ (jl.ip . 


(Step 1) The initial-correction. According to (16.7p . we take 4’o(3;) = x + a for some constant 
a, to be specified below. From ()6.5p and ()4.26p . it is not difficult to prove that 


( 6 . 8 ) 


TW+ln 4M'^21n^ + ln 

'^(x + 1) $o(a^ + l) 

4 + 


x + i 


x + a 
X + 1 + a 


+ o[^ 


Solve the equation —l/4 + o = 0, we get a = 1/4. By Mortici-transformation, we obtain 

( 6 . 9 ) 4 >o(x) = x + j, CFo(f(x)) = ^-Y- 

4 X + 4 
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As we need to use Mortici-transformation in each correction-process, so will not mention it for 
the sake of simplicity. 


(Step 2) The first-correction. Let us expand the following function into a power series in 
terms of 1 /x: 


( 6 . 10 ) 


In 


fix) 


d>o(x) + f 


/(i+i) 


= 21n 


X 


■1 . x + \ + ^ 


X + 


■ -|- In ■ 


_ 1 _ 

' x+1 


We solve the equation —1/16 + 2 ko = 0, and obtain kq = 1/32 / 0. Hence we take the 
first-correction MCi(a:) to be Type-I, i.e. 


( 6 . 11 ) 

Since 


MCi(x) 


X -|- Aq 


In ■ 


f (x) ^o(^) + x-l-Ao _ 128 32° _|_ Q 


f{x + l) ^o(a; + 1) + 
we enforce ^ = 0, and deduce Aq = |. Thus, 


x^ 


( 6 . 12 ) 


MCi(x) = 


32 


X + 


1 ’ 


CF.ifix)) = 


x + i + ^ 


(Step 3) The second-correction to the sixth-correction. Now we take MC 2 (x) to be 
Type-I, and let 


(6.13) 


MC2(x) 


^0 

^ + ^0 + 


By using (j4.26l) . we have 


(6.14) 


f{x) d'o(x)-bMC 2 (x) 

"" fix+ 1)^ ^ + 1) + MC2 (x + 1 ) 

512-1 , 5(-27 + 176ki + 64kiAi) . ^ 
x^ 2048x6 ^ \x^ 


Solve the equations 
(6.15) 


_9_ill — 0 

512 8 


—27 + 176ki + 64kiAi = 0, 
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we attain 


(6.16) 


Ki = 


64’ 


Ai — 


We take the /cth-correction MCfc(x) to be Type-I, then repeat the above approach like the 
second-correction, and solve successively the coefficients Hj and Xj (2 < j < 6) as follows: 


(6.17) 

(6.18) 


25 



49 


As 


1 81 1 
i’ ^4-^, A4--; 


K5 = 


121 

'm’ 



K6 = 


169 


Ae 


1 

4’ 


K7 = 


225 


A? — 


1 

4' 


From these results, it is not difficult to guess that 


(6.19) 




(2m -|- 1)^ 
64 


Am — 


1 

4' 


Further, we apply (j4.26p to check that the above conjecture holds true for some larger m. In 
this way, we hnally test that the fastest possible formula should be (16.3h . □ 


6.2 The continued fraction formulas involving the volume of the unit ball 

Let Xln be defined by dSH). The main purpose of this subsection is to present the following two 
theorems. 


Theorem 3. Let n > 1 be a positive integer. Then 


( 6 . 20 ) 






2(n + l) 


2n + l + K-=o 


/ (2m+l)n ■ 

2(2n+l) J 


Proof. . It follows from (j5.1l) and the recurrence relation r(3: + 1) = xT{x) that 

r(| + i)r(j + j) „ + ir'^(| + i) 

'■ ' r2(f + i) 2 r^d + i)' 

Replacing a; by ^ in (16.3|) . then after simplification, we get easily the desired assertion. 
Theorem 4. Let re G N, then 


□ 


( 6 . 22 ) 


1 

Hr). 


20F’ 


'2re + l+ K 


(2m-h 1)2 


'=0 \ 2(2re + 1) 


Proof. From (|5.ip . we have 
(6.23) 


Hn-l _ 1 r(| + 1) 
Hn V^F(^-|-^) 


Replacing x by ^ in (|6.3p . then taking reciprocals of both sides, finally substituting it into the 
above formula, this will complete the proof of Theorem 4. □ 
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Remark 9. Condition (|l.ll) is not an essential restriction. Actually, we can extend our method 
to any negative integer v. For example, by taking reciprocals of both sides in (|6.3p . we have 


(6.24) 


r^(g; + l) 

r2(x + i) 


1 1 oo 

= X + - H -— K 

4 X + 1+ m=l 


{2m+l)2 

64 

x + l 


Rex > 0. 


In this case, we take v = —1. It should be remarked that we can discover the above formula 
directly by using an approach similarly to Lemma 5. 

Remark 10. To the best of our knowledge, formula (16.ip and (j6..Sh were possibly neglected by 
many mathematicians for about more than twenty years, until 2013 I. Gavrea and M. Ivan 
mentioned it in their paper m- 


7 A continued fraction formula of Ramanujan 


The following lemma is Entry 39 in Berndt [8, p. 159], which is one of three principal formulas 
involving gamma functions given by Ramanujan. It is very difficult for us to imagine how 
Ramanujan discovered those beautiful continued fraction formulas. Maybe our method provides 
a theoretical basis. 

Lemma 6. Let I and n denote arbitrary complex numbers. Suppose that x is complex with 
Re X > 0 or that either n or I is an odd integer. Then 

^ ^ _ r (4(x + / + n + 1)) T (4(a; - / + n + 1)) r (4(x + / - n + 1)) T (4(x - / - n + 1)) 

T (l(x + / + n + 3)) T (|(x - / + n + 3)) r (i(x + / - n + 3)) T (l(x - / - n + 3)) 

8 1^-n^ 1^ -P 3^-n^ 3^ - P 

(x^ — P + n^ — l)/2 + 1 + x^ — 1 + 1 + x^ — 1 H- 

By replacing x by 4x, and taking {l,n) = (0,0), (/,n) = (1/4,1/2), {l,n) = (1/3,1/2), 
{I, n) = (1/8,1/2), respectively, the authors have checked that Lemma 6 is not optimal continued 
fraction expansion. Now, by employing these test, we may refine it in a uniform expression as 
follows. 


Theorem 5. Under the same conditions of Lemma 6, we have 


8 (1^ -n^)(l^-^^) {3^ - n^){3^ - P) 

^(x^ — P — + 1) — x^ — P — + (3^ + 1^ — 1) — x'^ — P — + (5^ + 3^ 

8 /—((2m — 1)^ — n^) ((2m — 1)^ —\ 

^(x^ — P — + 1) + m=i I x‘^ — P — + 8m^ + 1 j 


1 ) -••• 


Proof. We follow the method of Entry 25 in Berndt [8, p. 141]. Eirst, we rewrite Lemma 6 in 
the form 


8 1 G r)2 

| + ^(x2 + /2_n2-l) = x2-l+.^ 


1 


+ X 


12-/2 
2 - 1 


32 -n 2 


+ 


+ X 


32-/2 
2 - 1 


+ ■ 
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or 


1 1 12 - n 2 12-/2 32 _ ^2 32 _ i 2 

8/P+i(x2 + /2_n2-l) “ x2-l + 1 + x2-l + 1 + x2 - 1 + • • • ■ 

Secondly, by Entry 14 of Berndt [8, p. 121] (an infinity form see [8, p. 157]), we have 


8/P + i(x2 + /2 -n2 - 1) 

1 (l^-n^)(l^-/^) ( 3 ^-n^)( 3 ^-/^) 

a;2 — n2 — + (32 + 12 — 1) — — P — + (52 + 32 — 1) - 

((2m — 1)^ — n^) ((2m — 1)2 — /2) 

— x2 — /2 -72^-1- ((2m + 1)2 + (2m — 1)2 — 1) — .. . ' 


Note that (2m + 1)2 + (2m — 1)2 — 1 = Srn? + 1. Now take the reciprocal of both sides 
above and then solve for P, which again involves taking reciprocals. This will finish the proof 
of Theorem 5. □ 

The following theorem is the fastest possible form for Entry 26 in Berndt [8, p. 145]. 

Theorem 6. Suppose that either n is an odd integer and x is any complex number or that n is 
an arbitrary complex number and Rex > 0. Then 


r2 (4(x + n + 1)) (|(x - n + 1)) 
r2 (i(x + n + 3)) r2 (i(x - n + 3)) 

8 ^ —(2m — 1)^ ((2m — 1)2 — n^) 

^(x2 — n2 + 1) + m=i I x2 — n2 + 8m‘^ + 1 


Proof. Set / = 0 in Theorem 5, the desired equality follows at once. □ 

Similarly, we give another form of the Corollary in Berndt [8, p. 146]. 

Corollary 2. //Rex > 0, then 

r-(l(x+l)) 8 S- ( -(2m-1)* \ 

' ' rqiii + s)) i(i' + i) + miiVi^ + 8ni2 + iy 

Proof. We set n = 0 in Theorem 6, this completes the proof of the corollary readily. □ 


8 Some new conjectural continued fraction formulas 

In this section, we will give three examples to illustrate how to guess their fastest possible 
continued fraction expansions. Eor the recent results involving these functions, see Mortici, 
Cristea and Lu Cao and Wang and Chen m- 
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8 1 For 

o.± roi p3(3.+i) 


r3(x+i) 

In this subsection, we will use the function as an example to explain how to guess its 

fastest possible continued fraction expansion, which consists of the following steps. 


(1). Define 



Find the structure of CFk{f{x)) or MCk(x) by Mortici transformation and (I4.26p . We may 


determine that CF^ (fix)) has the form of Type-II, and its MC-point ui equals to 1/6. Here we 
omit the details for finding those coefficients in CF^ (/(x)), since the proof is very similar to 
that of Sec. 5 or Subsection 8.3 below. 


(2). We denote CF^ (/(x)) in the simplified form like (|2.6p : 


( 8 . 1 ) 




where A_i = 


(3). We write two sequences (Km)m>o and {Xm)m>o in the canonical form, then extract their 
common factors, respectively. For example, one may use Mathematica command “Factorinteger” 
to do that. In this way, we denote these two sequences in the form 



(4). Now we will look for the general terms of the sequences (Km)m>o and {Xm)m>o- We try 
to decompose them into some more simpler “partial sequences”. 

(4-1). We observe easily that the sequence {am)m>o = (1)3,5,...) has the general term 
am = 2m -|- 1. While, for the sequence (6m)m>o = (1 • 3,3 • 5, 5 • 7,...), its general term is 
bm = (2m -I- l)(2m -|- 3). 

(4-2). Let us consider the sequence (am)m>o = (1)2,4,5,7,8,10,11,...), which is the se¬ 
quence generated by deleting the sequence {3k)k>i from the positive integer sequence 
We can check that the sequence (am)m>o satisfies the following difference equation 



1, if m is an odd, 

2, if m is an even. 
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with the initial value uq = 1. Hence, we deduce that the general term equals to 


m 


= m + 1 + [yj. 


(4-3). Similarly, the sequence {/3m)m>o = (1,5,7,11,13,17,19,...) satisfies the following 
difference equation 


/3m Pm—1 — 


4, if m is an odd, 
2, if m is an even. 


with the initial condition /3o = 1, and we get 


f^m — ^Tn -|- 1 -|- 




(4-4). The sequence (Cm)m>o = (5^7, 3307,17167,5-31-353,5-13-2063, 5-7-19-419, 516847, 7- 
13 - 9697,...) is most difficult. Consider a new sequence (um)m>o to be defined by 

(8.2) Um'■= im mod (2m-|-l)(2m-|-3). 

By using Mathematica command “mod”, we can verify 

(8.3) [uq, ui, U2, - - -) =(58,220,490,868,1354,1948,2650, 3460,...) 

=2(29, no, 245,434, 677,974,1325,1730,...) 

:=2('Um)m>0- 

We may check that the sequence (um)m>o satisfies the following difference equation 


(8.4) Um 2Um—1 + Vra—2 — 108 

with the initial conditions vq = 29 and vi = 110. Solve this difference equation of order 2, we 
can deduce that the general term equals to 

(8.5) Um = 27(m + l)2 + 2. 

Now we rewrite the general term in the form 

( 8 . 6 ) Cm = 2(2m-M)(2m -h 3)um + w^m- 

We may check that {'Wo,'Wi,W 2 , ■ ■ ■) = (1,7,17,31,49,71,97,127,...). Quite similarly to the 
previous sequence (um)m> 0 ) {wm)m>o also satisfies a difference equation as follows 

(8.7) Wm + Win—2 = 4: 
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with the initial conditions tco = 1 and wi = 7. In this way, we get 


( 8 . 8 ) 


Wm = 2(m + 1)^ - 1. 


Substituting (8.5) and (8.8) into (8.6), we discover 
(8.9) im = 2(2m + l)(2m + 3) (27(m + 1)^ + 2) + 2(m + 1)^ - 1. 


Combining the above results and after some simplification, we conjecture that the general terms 
should be 


( 8 . 10 ) 

( 8 . 11 ) 




2 (m + l+[fj)^ (^3m + l + ^^^) 

7^ (2m +1)2 ’ 

A. (2(27(m + 1)^ + 2) + 

108 \ (2m + l)(2m + 3)/ 


(m > 0) 


(m > —1). 


Note that we used the fact that the last formula also holds true for m = — 1. 


(5). Define two sequences (fi;m)m>o and (Am)m>-i by (|8.10p and (18.lip , respectively. By 
making use of (14.2611 . we check that the above conjectures are still true for some “larger” m. 


(6). Further simplification for the general term Km- Actually, we have 


(8.12) (^m + l + lj\)(^3m + l + ((2m + 1)^ - (i)^^ 

(3m + l)(3m + 2)) 

“ 2 ’ 

which may be proved easily according to m is an odd and an even, respectively. Hence 


(8.13) 


— 


1 (3m +l)3(3m + 2)3 

2916 (2m + 1)2 


(m > 0). 


Finally, we propose the following reasonable conjecture. 


Open Problem 1. Let two sequences (Km)m>o and (Am)m>-i be define by (I8.13P and (18.lip , 
respectively. Let real x > —1/6, then we have 


(8.14) 


r3(x + i) 


3^ 


F3(x + 1) (x + 1)2 + A_1 + 






(a:+i)^+A2H-- 


(3;+4)^+A3 + 
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Remark 11. Open Problem 1 means that if there exists a fastest possible continued fraction 

r3(a;+i) 

expansion for the function then it must be the continued fraction expression of the 

right side in (j8.14l) . 


Replacing x by x — 1/6 and then after some simplification, we get the following equivalent 
forms of Open Problem 1. 


Open Problem Ih Let real x > 0, then 


(8.15) 


r3(x + |) 


, ^ , K 
® + 108 + *^=1 


(3n-2)3(3n-l)3 

2916(2n-l)^ 


^ + TM (2(27n2 + 2) + TnrzfYTo^ 




(2n-l)(2n+l) 


8.2 For 


r3(x+|) 


r3(3;+l) 

The main purpose of this subsection is to conjecture the fastest possible continued fraction 
expansion for the function /(x), which is defined by 

> r’(x + |) 

iWT) ■ 

We follow the same method described in last subsection. By testing, we observe that CP/ (fix)) 
has the form of Type-I, and its MC-point u is 1/3. Some computation data are listed as follows: 

1 


(8.16) 


where 


CF, {fix)) = 


x + i + K'^-ir 


\x+| 


(8.17) Ko = 


27’ 


(k1,K2,K3, 


1 /23 43 5^ 7^ 8^ 10^ ll^ 13^ 143 16^ 17^ 

M vT’y’ 


Similarly to the sequence {am)m>o in last subsection, it is not difficult to verify that the general 
term of the sequence {Xm)m>i should be 

18 181 A - 1 ("^ + 1 + LfJ)^ 

(8.18) ^^“ 54 2LfJ+l ■ 

Open Problem 2. For all real x > — |, we have 

r=^(x + |) 1 


(8.19) 


r3(x + i) 


X+ 3 + 


27 


31+3 + 


1 w 


1 4^ 

1 5^ 

1 , 

-7J +- 

3 :+^ + 
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Replace x by x — 1/3, we have the following equivalent forms of Open Problem 2. 


Open Problem 2 '. Let kq = and the sequence {Xm)m>i be defined as (IS.lSp . Let x > 0, 
then 


( 8 . 20 ) 


r^(x + ^) _ 1 /Km\ 

r3(x + |) X + m=0 \ X ) 


Since the partial coefficients of the continued fraction of the right side in (18.201) are all positive, 
we can prove the following consequence easily. 


Corollary 3. Let x > 0. Assume that Open Problem 2' is true, then for all non-negative integer 
k 


( 8 . 21 ) 


1 ^ r^(x + 3) ^ K 

X + j=o \ X ' r3(x + |) X + j=0 V X 


Remark 12. The authors have checked that Corollary 3 is true for k < 10. 


8.3 


JPqj. r(a:+??)r(3:+l-??) 

r2(3;+l) 


Let r; be a real number with 0 < rj < 1. In this subsection, we will discuss the continued fraction 
approximation for the ratio of the gamma functions 


( 8 . 22 ) 


Gr^ix) 


r(x + r 7 )r(x + 1 — ? 7 ) 
r2(x + i) 


It follows from (16.61) that 


(8.23) lim xGAx) = 1. 

X^OO ' 

Now let us begin to look for CFk{Gn{x)). 


(Step 1) The initial-correction. Note that 1 / = 1 in (j4.23l) . It follows readily from the 
recurrence formula r(x + 1 ) = xr(x) that 

Gn{x) _ (x + 1 )^ 

G^(x + 1) (x + r/)(x + 1 - 77 ) ’ 


Now we apply Mortici-transformation to determine <I>o(x) = x + cq. By making use of 
Mathematica software, one has 


(8.24) 


(x + r 7 j(x + 1 — 77 ) X + 1 + Co x^ \x J 


Solve Co — 77 + 772 = 0, we obtain co = 77 — 772 . 
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(Step 2) The first-correction. Let 


(8.25) 


MCi(x) = 


Ko 


X Xq 

similarly to the initial-correction, we also have 

(x-I-1)^ X-I-Co-I-MCi(x) 


(8.26) In 


-I- In ■ 


{x +'n){x + I - r]) X1Co-h MCi(x-h 1) 

2ko — + 2 ??^ — rj^ — 3^0 — 3koAo — S^oh + + 3ko?7^ — ‘irf — -q^ + 


x^ 


+ ^ + o 


X'^ 


1 


x^ 


where 


Solve 

(8.27) 

we get 

(8.28) 


g{x) =4ko - 2^0 + SkqAo -|- 4koAo -|- 6/tor/ -|- 4/toAor/ - 3r/^ - 2/tor/^ 

— 4/toAor/^ -|- 4r/^ — 8/tor/^ + 2//“^ -|- 4/tor/'^ — 2r/^ — 4r/® -|- 4r/'^ — q^. 

2/to — r/^ -I- 27/^ — q'^ = 0 

< —3/to — 3/toAo — 3/tor/ -|- 2r/^ -|- 3/tor/^ — 3r/^ — r/^ -|- 3r/^ — r/® = 0 

, ^(a;) = 0, 


(-1 r /) 2 r /2 

^0 — -- 


Ao 


1 — q + q'^ 
3 


(Step 3) The second-correction to the sixth-correction. Similarly to the first-correction, 
we use Mathematica software to find that the second-correction to the sixth-correction are the 
form of Type-I, and then solve all coefficients in these correction functions. It should be remarked 
that there is a parametric q, the related computations will become very huge and complex. So 
we need to manipulate Mathematica command “'Simplify ” . Here we list the final computing 
results as follows: 


Kl 

K2 

K3 

K4 

K5 


(-2 - q + r/^)^ 
36 ’ 

(—6 — q + r/^)^ 
100 ’ 

(—12 — q + r/^)^ 
196 

(—20 — q + q^Y 
324 

(-30 - q + r/^)^ 
484 


Ai — 
A2 = 


4 — 7/ -|- r/^ 


As — 



15 ’ 

9- 

7/ -1- 7/2 


35 ’ 

16 

— q + q^ 


63 

25 

— q + q^ 


99 

36 

-q + q^ 


143 
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Similarly to Open Problem 1, by careful data analysis and further checking, we may propose 
the following conjecture. 


Open Problem 3. For x > 0, then 


(8.29) 


r(x + r])r{x + 1 — 7]) 

r2(x + i) 


X + T]-T]^ + ( 




where kq 


i-vWf 

2 


and 


(8.30) 

(8.31) 


Krn 


( —m(m + 1) — 7/+ 7/^)^ (tTI + 7/)^(777 — ?7 + 1)2 
4(2777 + 1)2 - 4(2777 + 1)2 

(1 + 777)2 

-r, m> U. 

(2777 + 1) (2777 + 3)’ 


777 > 1 


Remark 13. If we take 7/ = ^, then the above conjecture implies (16.31) (i.e. the generalized Lord 
Brouncker’s continued fraction formula). Here we note that 


(777 + 77)2(777 — 77 + 1)^ (777 + i)^ (2777 + 1)2 

4(2777 + 1)2 “ 4(2777 + 1)2 “ 2^ ’ 

(1 + 777)^ - 77 + 77^ _ (1 + 777)2 - (i)2 ^ 1 _|_ 3^ 

(2777 + 1) (2777 + 3) (2777 + 1) (2777 + 3) (2777 + 1) (2777 + 3) 4 


9 Conclusions 

In this paper, we present a systematical way to construct a best possible finite and infinite 
continued fraction approximations for a class of functions. In particular, the method described 
in Sec. 4 is suitable for the ratio of the gamma functions, e.g. many examples can be found 
in the nice survey papers Qi [28] and Qi and Luo [29|. As our method is constructive, so 
all involving computations may be manipulated by a suitable symbolic computation software, 
e.g. Mathematica. In some sense, the main advantage of our method is that such formal continued 
fraction approximation of order k is the fastest possible when x tends to infinity. Concerning 
applications in approximation theory, numerical computation, our method represents a much 
better approximation formula than the power series approach (e.g. Taylor’s formula) for a kind 
of “good functions”. 

In addition, the multiple-correction method provides a useful tool for testing and guessing the 
continued fraction expansion involving a specified function. So our method should help advance 
the approximation theory, the theory of continued fraction, the generalized hypergeometric 
function, etc. Further, if we can obtain some new continued fraction expansions, probably these 
formulas could be used to study the irrationality, transcendence of the involved constants. 
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